Analogue Hawking Radiation in a dc-SQUID Array Transmission Line 
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We propose the use of a superconducting transmission line formed from an array of dc-SQUID 's 
for investigating analogue Hawking radiation. Biasing the array with a space-time varying flux 
modifies the propagation velocity of the transmission line, leading to an effective metric with an 
horizon. Being a fundamentally quantum mechanical device, this setup allows for investigations of 
quantum effects such as back-reaction and analogue space-time fluctuations on the Hawking process. 
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Introduction. — The possibility of observing Hawking 
radiation [l[ in a condensed matter system was first sug- 
gested by Unruh who uncovered the analogy between 
sound waves in a fluid and a scalar field in curved space- 
time 0]- In particular, the fluid equations of motion 
can formally be expressed in terms of an effective met- 
ric matching that of a gravitating spherical, non-rotating 
massive body in Painleve-GuUstrand coordinates Q 



v{rf] df + 2v{r)drdt + dr^ + r'^dif, (1) 



where Cg is the speed of sound and v{r) is the spa- 
tially varying velocity of the fluid. For a sound wave 
excitation in the fluid, with velocity Cs, the horizon oc- 
curs where f^(r) = and the excitation is incapable 
of surmounting the fluid flow. Since Unruh's original 
proposal. Hawking radiation analogues have been pro- 
posed using Bose-Einstein condensates N| , liquid Helium 
[5| , electromagnetic transmission lines j6l| , and flber-optic 
setups Estimated Hawking temperatures in these 
systems vary from a few nano-Kelvin to lO'^K respec- 
tively, far above temperatures predicted for astronomical 
black holes and thus usher in the possibility of experi- 
mental observation. Additionally, the understanding of 
the physics associated with laboratory system analogues 
may provide clues as to resolving unanswered questions 
associated with Hawking's original calculation such as 
the trans-Planckian problem [8] . 

In this letter, we propose using a metamaterial formed 
from an array of direct-current superconducting quan- 
tum interference devices (dc-SQUID's). Modulation of 
the propagation velocity, necessary for the formation of 
an horizon, is accomplished through application of an ex- 
ternal flux bias through the SQUID loops as indicated in 
Fig. [1^. Under appropriate conditions, this conflgura- 
tion provides the superconducting realization of Ref. , 
with the benefit of available fabrication methods. Indeed, 
arrays of SQUID's with parameters near those required 
to observe the Hawking effect have already been con- 
structed Furthermore, as a quantum device, the 
SQUID array goes beyond the capabilities of previously 
proposed systems, allowing the possibility to probe the 
effect on Hawking radiation of quantum fluctuations in 
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FIG. 1: a) Layout of the dc-SQUID transmission line. We 
assume each SQUID element is formed from identical tunnel 
junctions with critical current /c and capacitance Cj. b) Ef- 
fective lumped circuit model valid for frequencies below the 
plasma frequency and negligible SQUID self-inductance. 



the space-time metric. Thus, in principle, this setup en- 
ables the exploration of analogue quantum gravitational 
effects. 

Model. — We consider a coplanar transmission line 
composed of a centerline conductor formed by a long, 

3> 1, series array of dc-SQUID's indicated in Fig. [1^. 
For simplicity, we assume that all Josephson junctions 
(JJ) have identical critical current Ic and capacitance 
Cj values. For an individual dc-SQUID, with (pi and (/)2 
representing the gauge invariant phases across the JJ's, 
the equations of motion for 7± = {(j>i ± (f>2) /2 take the 
form 



1 d^7+ I 1 dj+ 
•jj'i, dt^ ojr dt 
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with plasma frequency Wp — (27r/c/Cj<i>o)^/^, character- 
istic frequency uJc = 2ttIcRn/^o, and normalized self- 
inductance /?l = 2ttLIc/^o- The parallel, normal current 
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resistance of the junction is denoted Rn, while $o = h/2e 
is the flux quantum and $ext is the external flux through 
the SQUID loop. If /3l < 1 then the SQUID dynamics 
can be approximated by a JJ with a flux-tunable critical 
current, = 2/c cos(7r$oxt/*I'o), the dynamics of which 
can be written 



1 



^+sm(7+) = - 



(3) 



where we have dropped the damping term, assuming the 
temperature is well below the superconducting critical 
temperature, and where the effective plasma frequency 
is given by uj^ — ^J2t:I^/{2Cj^o). We will assume the 
validity of this approximation and consider a flux-tunable 
array of Josephson Junctions (JJA). If we additionally 
restrict ourselves to frequencies well below the plasma 
frequency and currents below the critical current, then 
a JJ behaves as a passive, flux and current dependent 
inductance given by 



in(4,<fcxt) 
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(4) 



for the nth JJ in the array. The equivalent circuit is 
given in Fig. [TJd where we have labeled the length and 
capacitance to ground of each J J by a and Co, respec- 
tively. Using Kirchoff's laws, we can write the discrete 
equations of motion as 



ri+l 



dt 



In+l — In — 



dVn+ 

dt 



(5) 



From we see that by controlling the external flux 
bias, or by creating a varying current in the transmis- 
sion line, we are able to modify the inductance and thus 
propagation velocity inside the transmission line. Here, 
we focus on using the flux degree of freedom as our tun- 
able parameter. Creating a space-time varying current 
pulse, as in Ref. 0, can also be accomplished in our de- 
vice. However, our simplified model does not admit the 
correct dispersion relation to support the required sta- 
ble nonlinear solitonic localized pulses in the parameter 
region of interest. Charge solitons can however be pro- 
duced in the high impedance regime of our device 11]. 

Effective Geometry and Hawking Temperature. — By 
defining potentials An such that /„ — —CodAn/dt and 
Vn = An — An-1 [1], the equations of motion ^ can be 
combined to yield the discretized wave equation. 



dt dt 
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For wavelengths much longer than the dimensions of a 
single SQUID the dispersion relation becomes to lowest 
order in k: 
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where we have defined the velocity of propagation as 
c = a/y/LCo, which in practice is well below the vac- 
uum speed of light cq. In this limit, the wave equation 
approaches the continuum 
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By ignoring higher-order terms in Eq. ([7]), we effectively 
remove the discreteness of the array which, along with 
dispersion from JJ inertia terms, can play the role of 
Planck scale physics in our system 0, Sill]. For param- 
eter values considered below, the relevant short distance 
scale is c/ujp (> a). Requiring the propagation speed to 
vary in both space and time. 



^'{x — ut), 



(9) 



with fixed velocity u set by an external fiux bias pulse, 
the wave equation in the comoving frame becomes 
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A = 0, (10) 



where x and t now label the comoving coordinates. This 
wave equation can be re-expressed in terms of an effective 
space-time metric. 
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Comparing this metric with Eq. ([T]) , we see that our sys- 
tem contains a horizon located wherever v?' = c'^(x). In 
Fig. [2] we plot the effect of a step-like hyperbolic tangent 
fiux bias pulse of amplitude $ext = 0.2$o on a JJA with 
inductances given by Eq. (j4|) , where we have kept only the 
lowest term in the Ic/Ic expansion. Additionally, since 
$oxt can only increase the inductance, the flux-bias pulse 
velocity u must be below the unbiased transmission line 
propagation velocity c in order to establish a horizon. We 
do not consider Gaussian or similar pulse shapes as they 
generate both black hole and white hole horizons [13[ 
which complicates interpretation of the emission process. 

So far, we have focused on demonstrating a classical ef- 
fective background geometry with an event horizon. The 
next step is to quantize small perturbations in the po- 
tential field A about this background. The correct com- 
mutation relations between quantum field operators are 
required for conversion of vacuum fluctuations into pho- 
tons [IJ]. These relations have been verified in the sys- 



tems to which ours is analogous \B\ . The resulting Hawk- 
ing temperature is determined by the gradient of the JJA 
velocity at the horizon 



Th = 



2'Kkh 



dc{x) 



dx 



(12) 



The radiated power in the comoving frame coincides with 
the optimal rate for single-channel bosonic heat fiow in 
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axis of propagation (x) 

FIG. 2: Effect of a step-like flux pulse on the propagation 
velocity of a JJA as seen in the comoving frame. The pulse 
velocity was chosen to be u = 0.95c ($ext = 0). The black 
hole horizon occurs where c{x) = u. Arrow indicates the only 
permissible direction of travel across the horizon. 
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Eq. is universal for bosons since the channel- 

dependent group velocity and density of states cancel 
each other in one dimension [l^. For a detector at the 
end of the transmission line, the radiation emitted by 
an incoming bias pulse will be doppler shifted yielding 
higher power compared to Eq. However, the rate of 

emitted photons remains approximately unchanged. 

Model Validity. — For a single effective J J, the mag- 
nitude of quantum fluctuations in the phase variable 7+ 
depends on both the ratio of Josephson energy, Ej — 
$o-fc/27r, to charging energy, Ec = e'^/ACj, as well as 
on the impedance of the junction's electromagnetic en- 
vironment. These energy scales give a representation of 
the phase-charge uncertainty relation A-jAQ > e, and re- 
late the amplitude of quantum fluctuations between these 



variables [17|. When Ej/Ec > 1 and the impedance 



seen by the junction is less than the resistance quantum, 
Rq = /i/4e^ « 6.45kri, the phase operator behaves as a 
semiclassical quantity, i.e. the 7+ quantum fluctuations 
are small with respect to its average, and the JJ is in the 
superconducting state, allowing for a lumped inductor 
approximation. In the majority of experimental config- 
urations, a single JJ is connected to probe leads with 
impedance ^ 5051 and as such is in the low-impedance 
regime Z/Rq <C 1. In contrast, a JJA has an environ- 
ment that comprises not only the input and output ports, 
but also all the other JJ's in the array. In this case, we 
can define an effective impedance as seen by a single junc- 
tion to be Zj = Ze + Za where Ze is the environmental 
impedance of the leads and Za is the array impedance 




that, for freqtiencies below the plasma frequency, can be 
written as llSlI 



Za — Rq 



where the last equality explicitly shows the dependence 
on the external flux bias and single junction parameters. 
Thus, even for a small energy ratio Ec/Ej, the lumped 
inductor model applies only when Za/ Rq ^ 1- In Fig. [3] 
we show the dependence of array impedance Za on the 
external bias for fixed critical current Ic — 2 /iA and 
a range of experimentally valid capacitances to ground. 
As $cxt *i'o/2, high impedance causes large phase 
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FIG. 3: Ratio of array impedance Za to the resistance quan- 
tum Rq as a function of the external flux bias for a criti- 
cal current 7c = 2 jiA and example ground capacitor values: 
Co= 10"^*' (solid), 5 X 10~^''(dashed), 10""(dash-dot), and 
5 X 10~^*F (dotted). High and low impedance regions are de- 
fined above and below Za/Rq ~ 1 respectively. 



fluctuations, indicating a breakdown of our semiclassi- 
cal description; the array undergoes a quantum phase 
transition from superconducting to insulating Coulomb 
blockade behavior [l^. Note, the small J J parameter 
variability in actual arrays [l^, [l3| will prevent the di- 
vergence in Fig. [31 as well as cause some transmission line 
scattering in the low impedance superconducting state. 

The dependence of the Josephson energy Ej on exter- 
nal flux <I>oxt as described above allows for the systematic 
introduction of quantum fluctuations in our model. With 
the phase variable governing the circuit inductance, these 
fluctuations manifest themselves in the effective metric 
pT|) through the propagation velocity c. As the am- 
plitude of fluctuations increases, the metric becomes a 
quantum dynamical variable which must be included in 
the description of the Hawking process. Thus, conse- 
quences of back-reaction from the Hawking process as 
well as quantum dynamical space-time can be probed by 
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this configuration. Both processes, not included in the 
original Hawking derivation, represent analogue quantum 
gravitational effects present in our system [20|. 

Experimental Realization. — A possible realization of 
the JJA is shown in Fig. SI which consists of the JJA 
transmission line as well as an additional conducting line 
producing the space-time varying external flux bias $oxt- 
To provide a space-time changing velocity, the JJA is 
modulated by generating current pulses in the bias-line, 
the propagation velocity of which are assumed to be 
slightly below that of the unbiased JJA. The required 
bias pulse velocities u can be achieved by similarly em- 
ploying individual JJ's in series as the bias line. Ad- 
ditionally, a dc-external flux can be used to fine-tune 
the transmission line velocity closer to that of the bias- 
line, eliminating the need for large amplitude current bias 
pulses. Unavoidable current pulse dispersion in the bias- 
line, resulting in a decrease in Hawking temperature, can 
be minimized with appropriate choice of pulse shape and 
transmission line length. 



flux bias pulse 



ground plane 



SQUID array 



tunable phase-qubit 
detector 



-DODnaDDCl-A' » 1- 




FIG. 4: Possible transmission line and detector realization. 
Current pulses in the bias-line provide external flux necessary 
to modify the SQUID array propagation velocity. A phase- 
qubit at the end of the JJA functions as the photon detector. 



Unambiguous verification of the Hawking process will 
require frequency-tunable, single-shot photon detection 
at the end of the JJA opposite to that of the bias pulse 
origin. Although not presently available, microwave 
single-photon detectors based on superconducting qubits 
are under active investigation [2l|, . We will assume 
a phase-qubit as our model detector 



2l|. By repeat- 



edly injecting current pulses down the bias-line, the pre- 
dicted blackbody spectrum associated with the Hawk- 
ing process can be probed by tuning the qubit resonant 
frequency. Correlations across the horizon between the 
emitted photon pairs can be established through coin- 
cidence detection. We emphasize the essential need for 
correlation information in order to establish that a pho- 
ton is produced by the Hawking effect rather than some 
other ambient emission process, or spuriously generated 
via capacitive coupling to the bias-line. Unwanted direc- 
tional coupling can be minimized with proper engineering 
of the transmission line. 

To estimate the Hawking temperature we will assume 
parameters similar to those of Ref. with SQUID 's 
composed of tunnel junctions with Ic = 2 fiA and an 
upper bound achievable plasma frequency Wp = 27r x 
10^^ Hz. The capacitance to ground is assumed to be 



Co = 5 X 10-1^ F (dashed fine in Fig.^. Using a SQUID 
length a — 0.25 /im gives an unbiased transmission line 
velocity c ^ cq/IOO. Equation (fT^ gives the temperature 
as determined by the rate at which the JJA transmission 
line velocity varies that, in our case, is limited by the 
plasma frequency ujp. Assuming the maximum rate is 
an order of magnitude below cUp ($cxt — 0) /2tt, then the 
Hawking temperature is ~ 120 mK. This temperature 
can be a factor of ten larger than the ambient tempera- 
ture set by a dilution refrigerator and therefore should be 
visible above the background thermal spectrum. Using 
Eq. (fT3| and the sample pulse in Fig. ([2]) gives an initial 
Hawking temperature 120 mK, which decreases ~ 10% 
every 1000 JJA elements due to bias-line dispersion. Ap- 
plying the power expression (|13p yields an average emis- 
sion rate of one photon per pulse for ~ 4800 SQUID's. 
Of course, the transmission line can be made consider- 
ably shorter at the expense of an increase in number of 
pulse repetitions in order to accumulate sufficient photon 
counts to verify the Hawking radiation. The parameters 
and pulse shapes chosen here illustrate feasibility of this 
setup, but do not represent the only available configu- 
ration. These values can likely be improved upon and 
optimized in terms of both performance and fabrication 
of this proposal. 

Conclusion. — We have demonstrated that an array 
of dc-SQUID's in a coplanar transmission line, when bi- 
ased by a space-time dependent flux, creates an effective 
space-time metric with a horizon. As a quantum device, 
the superconducting transmission line allows for the pos- 
sibility of observing not only the Hawking effect, but also 
the effects of quantum fluctuations in an analogue grav- 
itational system. 
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